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SECTION  1 


INTRODUCnON 


1.1  BACKGROUND. 

Currentiv,  DELFIC  is  the  standard  code  for  prediction  of  radioactive  cloud  dispersion  and 
fallout.  It  is  sophisticated  in  the  treatment  of  processes  which  were  well  studied  at  the  time  of  its 
development  (the  1960s),  such  as  radioactive  d^y.  Its  treatment  of  hydrodynamics  and 
microphysics  is  more  crude,  however.  Research  in  hydrodynamics  and  microphysics  in  the  DNA 
and  other  communities  since  DELFIC  was  developed  has  advanced  the  state-of-the-art 
considerably,  particularly  m  the  field  of  turbulence.  It  was  therefore  recognized  that  first- 
principles  hydrodynamics  codes  developed  under  DNA  sponsorship,  such  as  MAZ  (Hatfield  et  al. 
1991),  could  be  upgraded  to  giv  more  accurate  predictions  of  radioactive  nuclear  clouds.  In 
particular,  microphysics  models  for  condc  isation  of  radioactive  species  and  agglomeration  of 
radioactive  particles  were  recognized  as  areas  where  improvement  was  needed. 

The  modeling  of  turbulence  in  the  simulation  of  nuclear  airblast  and  cloud  rise  has 
advanced  considerably  under  DNA  sponsorship  (see  Hassig  et  al.  1991,  1992;  Hookham  et  al. 
1991,  e.g  ).  Further,  the  phenomenon  of  turbulent  unmixedness  (Walitt  et  al.  1991)  has  been 
shown  to  be  potentially  important  when  computing  source  terms  for  non-linear  processes  such  as 
chemical  reactions.  The  conventional  approach  of  using  cell-averaged  values  of  turbulently 
fluctuating  quantities  to  compute  these  source  terms  can  give  erroneous  results.  Turbulent 
fluctuations  in  such  quantities  as  radioactive  vapor  concentration,  particle  velocity,  and 
temperature  could  have  a  significant  efiect  upon  microphysics  processes  such  as  condensation, 
agglomeration,  and  ablation.  The  efiect  of  turbulence  and  turbulent  unmixedness  upon 
microphysics  processes  was  therefore  identified  as  an  area  which  warranted  further  investigation. 

1.2  OBJECTIVES. 

The  primary  objectives  of  this  study  are  listed  below. 

•  Determine  the  sensitivity  of  the  distribution  of  radioactivity  within  a  nuclear  cloud  to 
agglomeration  and  condensation  models  and  their  parameters. 

•  Determine  which  physical  processes  are  important  to  model  and  which  are  not. 

•  Perform  sensitivity  studies  on  parameters  which  are  not  well  characterized. 

•  Develop  improved  condensation  and  agglomeration  models  for  radioactive  species. 

•  Investigate  the  effects  of  turbulence  upon  condensation  and  agglomeration. 
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O  APPROACH. 


The  sensitivity  of  the  caloilated  distribution  (spatial  and  temporal)  of  radioactivity  in  a 
nuclear  cloud  to  agglomeration  and  condensation  models  was  studied  by  performing  parametric 
variations  on  a  test  case,  described  below. 

Test  Case 

•  300  kt  at  50  sft 

•  Production  zoning  (10,000  zones  max.) 

•  6  particle  size  groups;  mean  diameter  ==  4.S,  32,  1 12,  and  500  ^  and  0.2  and  0.80  cm  for 
groups  1-6 

•  One  composite  radioactive  species  based  on  curve  T  McGahan  (1990) 

•  Current  standard  nuclear  dust  cloud  physical  sub-models 

•  Moscow  July  atmosphere  (relatively  humid) 

The  test  case  was  re-run  as  changes  in  the  agglomeration  and  condensation  models  were 
made.  The  calculations  were  performed  with  the  MAZ  hydrocode  (Hatfield  et  al.  1991). 
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SECTION  2 


RADIOACTIVE  DUST  CLOUD  MICROPHYSICS  MODELS 

The  microphysics  model  improvements  wilt  be  separated  into  "basic"  model  improvements 
and  turbulent  unmbcedness  model  improvements.  The  e^ect  of  the  model  changes  will  be 
illustrated  by  comparing  results  &om  the  test  case  as  each  change  was  made.  The  matrix  of  test 
cases  and  model  changes  is  shown  in  Table  1 . 

2.1  BASIC  MODEL  IMPROVEMENTS. 

The  "baseline"  model  for  condensation  of  radioactive  vapor  and  subsequent  agglomeration 
of ''adioactive  particles  was  adopted  to  give  a  preliminary  estimate  of  the  distribution  of 
radioactivity  in  a  nuclear  cloud.  The  model  considered  agglomeration  between  pure  radioactive 
particles  and  non-radioactive  dust  only,  and  neglected  radioactive  dust-dust  agglomeration  and 
pure-pure  agglomeration.  In  an  attempt  to  study  the  effect  of  agglomeration,  the  radioactive 
vapors  were  forced  to  condense  preferentially  upon  the  pure  radioactive  particles  relative  to  dust 
particles  in  a  90/10  ratio  (if  this  were  not  so,  most  of  the  vapor  would  condense  upon  the  dust, 
which  was  not  allowed  to  agglomerate).  A  1  |xm  size  was  assumed  for  the  pure  condensed 
radioactive  particles  (conrads).  Case  TMOl  used  these  baseline  models. 

The  first  model  improvement  was  to  allow  radioactive  dust-dust  agglomeration  and  to  not 
force  the  radioactive  vapors  to  condense  preferentially  upon  the  pure  particles.  Agglomeration 
was  allowed  to  occur  as  a  result  of  relative  mean  velocities  between  particles  of  different  sizes, 
and  thus  since  all  particles  of  a  given  size  have  the  same  velocity,  particles  of  the  same  size  could 
not  agglomerate.  The  agglomeration  model  for  this  "differential  scavenging"  process  is  described 
in  Figure  2-1. 

Calculation  TM02A  included  dust-dust  agglomeration  and  removed  the  restriction  on 
where  radioactive  vapors  could  conden.'!';.  Calculation  TM02B  investigated  the  effect  of  lowering 
the  particle  collision  efficiency  to  1  fi'om  the  previous  value  of  2  (2  was  used  previously  to 
account  for  the  effects  of  turbulence  upon  agglomeration  in  a  very  crude  wayi  Calculation 
TM02C  then  investigated  the  effect  of  changing  the  "sticking  efficiency"  from  the  nominal  value 
of  1  (an  upper  bound,  since  presumably  all  particles  that  collide  will  not  stick)  to  1  if  the  particle 
surface  temperature  was  greater  than  the  melting  point  of  dust,  and  0  if  it  was  not.  The  latter 
condition  is  probably  closer  to  a  lower  bound  on  sticking  efficiency.  TM03  then  added  a  simple 
turbulent  agglomeration  model,  described  in  Figure  2-1,  which  accounted  only  for  agglomeration 
of  similar-size  particles  which  were  smaller  than  the  small  turbulent  eddies  in  the  flow  (there  is  no 
such  simple  model  for  larger  particles  which  do  not  follow  the  flow,  although  we  can  addrer.s  this 
issue  with  our  unmixedness  model). 

The  distribution  of  radioactivity  within  a  dust  cloud  can  be  effectively  summarized  by  the 
cumulative  activity  particle  size  distribution  (APSD)for  the  cloud,  which  is  shown  in  Figure  2-2  for 
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Figure  2-1.  Agglomeration  model  description 
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Figure  2-2.  Cumulative  activity  particle  size  distribution  comparison 


the  e&i^  listed  in  Table  1  (all  curves  are  for  1  minute,  except  TM03A,  which  is  shown  at  1  and  10 
minutes).  Case  TMOl  shows  a  large  amount  of  pure  radioactive  particles,  which  agglomerate 
relatively  rapidly  due  to  their  small  size.  (Agglomeration  rate  due  to  particle  relative  mean 
velocity  is  inversely  proportional  to  particle  radius,  neglecting  the  effect  of  lower  collision 
efficiencies  for  small  particles,  which  can  flow  around  larger  particles  instead  of  colliding.)  This 
calculation  is  shown  for  comparison  only,  as  it  does  not  represent  a  realistic 
condensation/agglomeration  model.  Proceeding  to  Case  TM02A,  allowing  dust«dust 
agglomeration  and  not  forcing  the  bulk  of  the  radioactivity  to  condense  on  the  conrad  (pure 
radioactive)  particles  resulted  in  a  reduction  in  conrad  concentration  from  ~80  to  -'30%  of  total 
activity.  The  amount  of  activity  on  the  larger  dust  particles  increased  also.  Lowering  the  collision 
efficiency  to  1  in  Case  TM02B  resulted  in  somewhat  less  agglomeration,  and  hence  less  activity 
on  the  larger  particles. 

Assuming  a  sticking  efficiency  of  1  for  hard  dust  particles  seemed  intuitively  to  he  an 
overestimate.  Thus,  the  sticking  efficiency  was  made  a  function  o.^  the  presence  of  physical  forces 
that  would  act  to  hold  the  particles  together.  Making  the  sticking  efficiency  a  function  of  whether 
liquid  is  present  on  particle  surfaces  (Case  TM02C)  resulted  in  a  substantial  decrease  in 
agglomeration  rate,  and  only  ~2%  of  the  activity  resides  on  larger  particle  sizes.  Allowing 
agglomeration  when  there  was  condensed  water  present  (Case  TM03A)  increased  agglomeration 
substantially,  particularly  if  the  calculation  was  continued  to  10  minutes.  The  effect  of  wet 
agglomeration  is  further  illustrated  in  Figure  2-3,  which  shows  vertical  distributions  of  cumulative 
attached  activity  and  agglomerated  activity  (pure  radioactive  particle  to  dust  agglomeration  only 
is  shown)  for  all  dust  groups  together  and  individually.  Wet  agglomeration  clearly  makes  a 
substantial  difference  in  the  amount  of  agglomeration,  at  least  for  the  relatively  humid  atmosphere 
used  for  these  calculations. 

Adding  the  simple  turbulent  agglomeration  model  in  Case  TM03  resulted  in  essentially  no 
change  in  the  APSD.  We  must  remember,  however,  that  this  model  only  considers  small  particles 
moving  with  the  turbulent  eddies,  and  neglects  entirely  larger  particles  moving  relative  to  the 
turbulent  eddies. 

In  addition  to  the  effect  of  agglomeration  on  the  radioactive  APSD,  adding  agglomeration 
of  non-radioactive  dust  caused  significant  changes  in  the  non-radioactive  dust  cloud.  Figure  2-4 
shows  a  comparison  of  Case  TMOl,  which  had  no  dust-dust  agglomeration,  to  Case  TM03,  which 
agglomerated  if  particles  were  melted,  and  TM03A,  which  agglomerated  if  particles  were  melted 
or  wet.  Case  TM03  gave  similar  results  to  Case  TMOl,  but  Case  TM03A,  while  having  a  similar 
amount  of  total  dust  lofred,  lofted  a  substantially  larger  amount  of  the  larger  size  particles  (group 
6,  e.g.).  (Note  that  in  the  lower  figures  the  numbered  contours  correspond  to  size  group  number, 
while  in  the  upper  figures  they  refer  to  concentration.) 

As  an  indication  of  how  the  APSD  varied  with  location  in  the  dust  cloud.  Figure  2-5 
shows  APSDs  computed  for  horizontal  "fly-throughs"  at  several  locations  for  Case  TM03A.  The 
fly-through  calculations  simulate  an  aircraft  flying  a  straight  horizontal  path  through  the  cloud  and 
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Total  Solid  Dust  Density  Contoufs  at  1  Minute 


collecting  mass.  It  is  evident  that  the  doud  st^  contains  many  more  large  particles  than  the 
cloud  as  a  whole,  while  the  top  of  the  cloud  contains  more  small  particles. 

2.2  TURBULENT  UNMIXEDNESS  CONDENSATION  MODEL. 

The  condensation  model  was  chosen  as  a  first  test  of  turbulent  unmixedness  modeling. 
The  model  is  described  in  some  detail  below. 

Consider  condensation  of  a  trace  quantity  of  radioactive  vapor  species,  with  mass  fi'action 
Cg,  on  solid  particles  with  number  density  .  In  order  to  evaluate  the  time  rate  of  change  of 
mean  vapor  concentration,  we  need  to  know  the  mean  condensation  rate.  Figure  2>6  expresses 
the  instantaneous  rate  as  a  function  of  local  gas  and  particle  properties.  The  relationship  is  non* 
linear  in  the  vapor  and  particle  concentrations,  temperature,  and  velocities.  Hence,  the  mean 
condensation  rate,  say  averaged  over  a  timestep,  cannot,  in  general,  be  evaluated  by  using  only 
mean  cell  properties  (see  Figure  2-6).  If  we  perform  a  Taylor  series  expansion  on  the  rate 
expression,  we  obtain  equation  2-4  in  Fig.  2-7,  indicating  a  second-order  concentration 
correlation.  (The  dots  in  the  equation  refer  to  additional  correlations  in  the  expansion.)  This 
correlation  can  be  positive  or  negative.  When  for  example,  at  interfaces  where  intermixing  of 
vapor  and  solid  particles  (condensation  sites)  occurs,  the  correlation  is  negative,  indicating  a 
condensation  rate  less  than  that  predicted  using  conventional  (mean-value)  models.  This 
phenomenon  is  referred  to  as  unmixedness  in  turbulent  combustion. 

To  proceed,  we  require  the  second-order  concentration  correlation.  Following  standard 
turbulent  closure  model  practices,  we  develop  a  continuity  equation  for  the  correlation  by 
mathematically  manipulating  the  individual  vapor  and  particle  continuity  equations  and  taking  an 
ensemble  average  of  the  result.  Sparing  the  details.  Figure  2-8  shows  the  final  equation,  in  which 
triple  correlations  have  been  replaced  by  modeled  diffusion  and  dissipation  terms.  These  modeled 
terms  are  analogous  to  similar  terms  in  the  air  closure  equations.  For  small  particles,  which 
follow  local  air  eddies,  we  set  the  empirical  coefficients  in  the  modeled  terms  to  the  air  values, 
which  have  been  calibrated  to  experiment.  For  larger  particles  which  slip  relative  to  the  gas 
phase,  the  coefficients  may  need  adjustment.  The  research  of  Rizk  and  Elghobashi  (1989),  which 
indicates  particle  size  dependency  on  turbulent  diffiision  velocities,  will  be  of  value  here. 
Ultimately,  the  validity  of  models  is  established  by  matching  calculations  with  controlled 
experiments. 

The  last  term  in  equation  2-S  represents  the  high-order  moments  of  the  condensation  rate 
in  the  vapor  continuity  equation.  Conventional  closure  methods  do  not  permit  its  evaluation. 
Instead,  we  apply  the  turbulent  chemistry  pdf  approach.  Figure  2-8  shows  2-D  slices  of  the  pdf 
for  vapor  and  particle  concentrations.  (The  pdf  is  actually  a  surface  in  c^-N  space.)  The  pdf  can 
be  approximated  with  a  finite  number  of  Dirac  fiinctions  with  as-yet  undetermined  strengths  e^. 
Constraints  include  the  values  of  mean  concentrations  and  double  correlation.  (In  an  explicit 
numerical  scheme,  we  take  values  fi-om  the  previous  timestep  as  known  quantities  in  order  to 
advance  these  quantities  to  the  next  cycle.)  The  locations  of  the  Dirac  functions  are  specified, 
discrete  events;  they  are  chosen  to  represent  the  physics  of  the  problem. 
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RADIOACTIVE  VAPOR  CONDENSATION 
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Figure  2-6.  Radioactive  vapor  condertsation  on  solid/liquid  particles. 
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s  IS  analogous  to  unmixedness  effects  on  gaseous  chemical  reactions  In  turbulent  flows. 
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Figure  2-8.  Method  of  solution:  radioactive  vapor  condensation. 
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Figure  2-8.  Method  of  solution;  radioactive  vapor  condensation  (Continued) 


Suitable  events  should  include  states  at  the  mean  {i.e.,  when  turbulence  levels  are 
lowl  as  well  as  states  where  one  or  both  concentrations  disappear  (intermittency). 
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Figure  2-8.  Method  of  solution:  radioactive  vapor  condensation  (Continued) 


For  this  appUeation,  we  choose  three  events  as  indicated  in  Figure  2-8.  One  state  is 
naturally  chosen  at  the  mean  concentrations  to  get  a  reasonable  state  description  when  turbulence 
levels  are  low.  Local  fluctuations  can  force  concentrations  of  one  or  both  species  to  zero 
(intermittency),  hence.  Dirac  functions  are  introduced  at  these  states.  Instantaneous 
concentrations  greater  than  the  mean  are  accommodated  with  Dirac  functions  placed  at  a  multiple 
of  the  mean  concentration.  Constraint  equations  2-6  through  2-8  allow  the  exact  computation  of 
Dirac  fiinction  strengths  for  3  states.  More  states  may  be  desirable,  however.  In  this  case,  there 
are  more  unknown  quantities  than  constraints,  so  minimization  of  an  "entropy"  function  "F",  can 
be  employed.  Note  that  for  the  test  calculation  described  below,  the  turbulent  correlation  CgN^ 
was  calculated  using  a  mixing-length  approximation  to  determine  if  turbulent  unmixedness  was 
important  to  the  condensation  process  before  the  full  second-order  closure  correlation  was 
calculated. 

Once  discrete  states  have  been  chosen  (in  each  computational  cell  based  on  local  physics), 
and  the  Dirac  function  strengths  are  determined  (see  equation  2-9,  which  follows  from  equations 
2-6  to  2-8  and  states  indicated  in  Figure  2-8),  any  high  order  moment  involving  vapor  and/or 
particle  concentrations  can  be  computed.  Equation  2-10  in  Figure  2-8  shows  an  example  for  a 
triple  correlation  which  is  needed  to  close  the  double  concentration  correlation,  equation  2-S. 

Figure  2-9  shows  a  comparison  of  Case  TM04,  which  used  the  turbulent  condensation 
model  described  above,  to  Case  TM03,  which  did  not  (the  calculations  were  otherwise  the  same). 
It  can  be  seen  that  adding  the  turbulent  unmixedness  model  had  only  a  small  effect  upon  the 
amount  of  condensation.  Examination  of  the  solution  revealed  that  although  the  unmixedness 
model  did  change  the  instantaneous  condensation  rate  substantially,  the  overall  rats  of 
condensation  was  limited  by  the  rate  of  cooling  of  a  given  computational  cell,  so  that  increasing 
the  instantaneous  rate  had  a  small  effect. 


SECTION  3 

CONCLUSIONS  AND  RECOMMENDATIONS 
3.1  CONCLUSIONS, 

The  conclusions  from  this  study  are  summarized  as  follows 

•  Adding  agglomeration  of  wet  particles  (mud)  produced  a  substantial  amount  of 
agglomeration  which  continued  to  about  10  minutes  This  phenomenon  had  the  single 
biggest  effect  on  agglomeration. 

«-  Basing  the  agglomeration  model  only  on  whether  the  particle  surface  was  melted  resulted 
in  a  relatively  small  amount  of  agglomeration.  If  a  sticking  efficiency  of  1  was  assumed 
even  when  a  particle  was  not  melted,  significant  agglomeration  occu.ted. 
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SECTION  3 


CONCLUSIONS  AND  RECOMMENDATIONS 

3.1  CONCLUSIONS. 

The  conclusions  from  this  study  are  summarized  as  follows. 

«  Adding  agglomeration  of  wet  particles  (mud)  produced  a  substantial  amount  of 
agglomeration  wluch  continued  to  about  10  minutes.  This  phenomenon  had  the  single 
biggest  effect  on  agglomeration. 

•  Basing  the  agglomeration  model  only  on  whether  the  particle  surface  was  melted  resulted 
in  a  relatively  small  amount  of  agglomeration.  If  a  sticldng  efficiency  of  1  was  assumed 
even  when  a  particle  was  not  melted,  significant  agglomeration  occurred. 

•  The  simple  Saffinan-Tumer  agglomeration  model  did  not  cause  a  substantial  amount  of 
agglomeration  in  addition  to  that  caused  by  relative  mean  particle  velocity  (differential 
scavenging). 

•  The  addition  of  a  turbulent  condensation  model  produced  only  a  small  change  in  the 
radioactive  cloud,  since  condensation  was  generally  limited  by  the  rate  of  cooling  in  the 
cloud,  rather  than  local  radioactive  vapor  and  particle  concentrations  fluctuations. 

•  Agglomeration  of  non-radioactive  dust  had  a  small  effect  on  the  total  dust  mass  in  the 
cloud,  but  increased  significantly  the  amount  of  larger  particles,  assuming  particle  sticking 
based  on  melted/not  melted  or  wet.  If  a  sticking  efficiency  of  1  is  assumed  for  all 
conditions,  total  mass  in  the  cloud  at  1  minute  drops  by  about  a  factor  of  three  and  the 
amount  of  larger  particles  increases  further. 

3.2  RECOMMENDATIONS. 

It  is  recommended  that  the  full  turbulent  unmixedness  model  for  agglomeration  be  tested. 
Although  the  simple  turbulent  agglomeration  model  did  not  show  an  effect,  that  does  not 
necessarily  mean  that  the  full  model  will  not,  since  the  simple  model  was  limited  in  scope  (the 
simple  model  only  considered  small  particles  moving  ^^dth  the  turbulent  eddies,  whereas  in  the  real 
case  there  win  be  substantial  numbers  of  larger  particles  which  do  not  follow  the  flow). 
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